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Abstract 

Self-consistent excited states of condensates are solutions of the Gross-Pitaevskii (GP) equation 
and have been amply discussed in the literature and related to experiments. By introducing a 
more general mean-field which includes the GP one as a special case, we find a new class of self- 
consistent excited states. In these states macroscopic numbers of bosons reside in different one- 
particle functions, i.e., the states are fragmented. Still, a single chemical potential is associated 
with the condensate. A numerical example is presented, illustrating that the energies of the new, 
fragmented, states are much lower than those of the GP excited states, and that they are stable to 
variations of the particle number and shape of the trap potential. 

PACS numbers: 03.75.Hh,03.65.Ge,03.75.Nt 
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Numerous properties of dilute Bose-Einstein condensates are well discussed by the Gross- 
Htaevs.„ ,GP) e,ua«„„ QQ. For .ev.ew. see. e.g. BQQ. K *e co,*,..e is ..pped. 

this equation, which is equivalent to the so called nonlinear Schrodinger equation, possesses 
a discrete spectrum of stationary states. Being a nonlinear equation which has to be solved 
self-consistently, the solutions of the GP equation are generally called self-consistent states. 



The properties of GP sel 
discussed, see, e.g. 



consistent excited states and their formation have been amply 
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known self-consistent excitations, such as vortex states 
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22, |2J,|2J|. The observation of solitons in trapped Bose gases 
provides a striking manifestation of nonlinear atom optics 

By introducing a more general mean-field which includes the GP one as a special case, 
we find another class of self-consistent excited states. Numerical examples show that the 
corresponding energies are substantially lower than those of the GP excited states. To avoid 
misunderstanding we mention that self-consistent excited states are physically distinct from 
collective (or particle-hole) excitations. The latter correspond to small oscillations around a 
given state and are described by mean-field /mear-response theories based on the Bogoliubov 



approximation 



We consider a system of N identical bosons. The corresponding GP equation is the 
mean-field equation of this system with a 5-function interaction potential W{fi — Tj) = 
Ao5(^j — rj), where is the position of the i-th boson and the nonlinear parameter Aq is 
related to the s-wave scattering length of the bosons |^. The GP ansatz further assumes 
the wave function \1'gp to be a product of identical spatial one-particle functions (orbitals) 
ip: '^Gp{^i:^2, ■ ■ ■ j^n) = f{ri)(f{r2)---f{^iy). With this wavefunction the total energy 
Egp =< ^gpI-^I^gp > is obtained as the expectation value of our Hamiltonian H and 



reads 



Egp = N{ ip*hipdf+ 



Ao(iV 



\ip\^df}, 



(1) 



where h{r) = T + V{r) is the unperturbed one-particle Hamiltonian consisting of the trap 
potential V{r) and the kinetic energy operator T. By minimizing the energy (1) one obtains 
the GP equation: 

{h{r) + Xo{N-l)\^{r)\'}^{r}= f^GPVi^. (2) 
Each self-consistent solution of this eigenvalue equation determines (p and the chemical 
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potential /icp- The ground state of the system, of course, corresponds to the solution with 
the lowest energy Egp- The other solutions describe the self-consistent excited states of the 
system. 

In the absence of the mutual interaction between the bosons, all bosons reside in a single 
spatial orbital in the ground state of the system. The GP ansatz for the wavefunction \1/gp 
is thus very appealing for the ground state and has indeed been very successful in explaining 
many observations. We have reason to assume, however, that not all relevant macroscopic 
excited states are describable by "^gp- The general wavefunctions \1/ of N non-interacting 
bosons is a product of orbitals which can all be different. Since the bosons are identical, this 
product must be symmetrized and several bosons may reside in the same spatial orbital. 
We may put rii bosons in orbital 0i, n2 in orbital 02 and so on. For transparency of 
presentation we restrict ourselves in the following to two orbitals 0i and 02 with particle 
occupations rii and ^2 and ni + n2 = N. The extension to more orbitals is straightforward. 
The wavefunction now reads 

^'(fi, . . . , Tat) = 50i(fi) ■ ■ ■ 0i(f„j02(r;^+i) ■ ■ ■ 02(r;^+„2), (3) 

where S is the symmetrizing operator. 

By definition the most general mean- field en ergy is E =< > where \I' is the above 

discussed wavefunction. It is evaluated to give 30| 

E = nihil + n2h22 + \o ^ / \(pi\^df + 

^^n2(n|-2) I |0^|4^^-.+ 2Aonin2y" \<Pi\'\<P2\'dr, (4) 

where ha = J (p* h(l)idf is the usual one-particle energy. We now minimize this energy with 
respect to the orbitals under the constraints that they are orthogonal and normalized, i.e., 

< 0l|02 >= 0,< 0i|0i > = < 02102 >= 1- (5) 

This leads to the following set of coupled equations for the optimal orbitals: 
{ h{r) + AoK - l)|0i(f)|2 + 2Aon2|02(r)|'} 0i(f) = 

= /ill 01 (^) +/il2 02(r) 

{ h{f) + AoK - l)|02(r)P + 2Aoni|0i(f)|2} 02(f) = 

= /i22 02(r) +/X210l(^- (6) 



These equations should not be confused with those for condensates made of two types of 
bosons. The fiij are the Lagrange parameters due to the above mentioned constraints; fin 
and /X22 are related to the normalization constraints and the introduction of /ii2 and /i2i 
enforces orthogonality of (pi and 02- We mention that nifii2 = n2fi2i- In general, /ii2 7^ 
unless 01 and 02 are of different spatial symmetry. 

By setting n2 = we readily see that the general mean-field wavefunction \1/ in (3) 
contains the GP wavefunctions "^gp as a special case. The energy E in (4) then reduces to 
the expression for Eqp in (1). Because of ^1/112 = n2/i2i, the parameter /ii2 vanishes and the 
first equation in (6) is then nothing but the GP equation (2). Contrary to the GP equation 
(2), the general mean-field equations (6) are - owing to the presence of /ii2 and fi2i - not 
eigenvalue equations. This is a generic property of the latter equations. Unless either one of 
the particle occupations vanishes, or rti = n2 = N/2, or 0i and 02 have different symmetry, 
or /ill = Ai22, the off-diagonal Lagrange parameters /ii2 and H21 cannot be removed from (6) 
by linear transformations. 

While we cannot attribute a physical observable to the off-diagonal Lagrange parameters 
/ii2 and /i2i we can do so for the diagonal ones. As can be seen from (1), the energy needed 
to remove a boson from the condensate without changing the orbital (p is, within GP mean- 
field, given by Egp{N) — Egp{N — 1) = ficp- Analogously, we can compute within the 
present general mean-field the energy needed to remove a boson from orbital 0i and that 
from orbital 02- Recalling that 0i and 02 are orthogonal, we readily find from (4) and (6) 
the identities /in = E{ni,n2) — E{ni — 1,^2) and /i22 = E{ni,n2) — E{ni,n2 — 1). Clearly, 
/ill and /i22 can be viewed as chemical potentials of the {0i}- and {02}-boson manifolds, 
respectively. 

Next, we can make use of the fact that the energy E and the orbitals 0i and 02 depend on 
the particle occupation ni and treat this occupation as a variational parameter, i.e., search 
for its optimal value which makes the energy stationary. Interestingly, one can show for 
macroscopic occupancies ni,n2 ^ 1 that at any extremum of the energy as a function of rii 
the two quantities /in and /i22 coincide j^^: 

/ill = /i22- (7) 

This is a relevant finding for understanding the concept of a condensate and the meaning 
of its self-consistent states. The system described by the wavefunction (3) consists of two 

4 



(or more) subsystems each possessing its own chemical potentiaL This may contradict the 
picture one usually has of a condensate. It is thus important to note that these, in general 
different, chemical potentials become identical at the optimal occupations restoring thereby 
the picture of a condensate. With this in mind, we define as self-consistent condensate states 
of the present mean-field only those solutions of (6) at which the chemical potentials /in 
and /i22 are identical. These correspond to the extrema of the energy E as a function of the 
boson occupation rii. 

Let us refer here to the problem of fragmentation often discussed in the literature, see 
.efs. a a„d .eferences tKe.ein. A condensate ,s f„ed ,f its .edueed one-body 

density matrix has two or more macroscopic eigenvalues (the total number of particles N 
is assumed large). Until now fragmentation has not been found for trapped condensates. 
Our present ansatz has the potential to describe fragmentation on the mean-field level. If 
the optimal boson occupations Ui and ^2 are both macroscopic for large N, fragmentation 
in the respective self-consistent state indeed takes place. Consequently, we call such states 
which are beyond reach of the GP equation (2), self- consistent fragmented states. 

We shall demonstrate in the following that the present mean-field supports self-consistent 
fragmented excited states. Moreover, these states can be at much lower energy than all the 
GP self-consistent excited states. As an example we choose a repulsive condensate (Aq > 0) 
in an one-dimensional double-well potential. Since equations (6) have very recently been 
derived, the numerical procedures to evaluate them are not yet sufficiently developed to 
allow for computations in three dimensions. We have reason to believe, however, that 
similar results will be obtained in more dimensions. 

In our example the potential 

V{x) = c + |(a;2 - ^/a^ + (A - 2xxo)^) (8) 

describes two wells separated by a barrier. The depths of these wells differ by a bias ~ A. 
The parameter c is chosen to set the bottom of the potential equal to zero. The coordinate x 
is dimensionless and all energies and Aq are now in units of the frequency u. The equations 
for the optimal orbitals have been evaluated numerically by using the DVR method 
Starting from an initial guess for (pi and 02, for instance, the solutions of the equations 
with shghtly different boson occupations, the equations are iterated until self-consistency is 
achieved. To accelerate the convergence of the calculations we have employed an energy- 
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shift technique as often done in Hartree-Fock electronic structure computations js^. The 
GP equation has been evaluated by solving the corresponding eigenvalue equation self- 
consistently using the DVR method. 

Having obtained the optimal orbitals for the present and GP mean-fields, we have com- 
puted the energies E and Ecp. Computations have been carried out systematically for 
the ground and the first few excited states using various values of the coupling constant 
(A = AqA^) (note that ^ 1) and of the parameters appearing in the potential V{x) in 
(8). For each set of these values the optimal orbitals and the mean-field energy E have been 
determined as a function of the boson occupation rii. A typical example is shown in fig.l. In 
this figure the mean- field energy per particle E/N is depicted as a function of the fractional 
occupation rii/N for A = 2.5. For comparison, the energies per particle Egp/N of the GP 
ground and first excited state are indicated as well. Two disjoint E/N curves are seen in 
the figure and we concentrate first on the one at lower energy. Starting the calculations at 
rii = N, the energies E and Egp, of course, coincide. We then continuously lowered the 
value of rii/N and obtained a smooth E/N curve. This curve exhibits two maxima and a 
minimum which is, of course, located between them. Following (7) there is only one value 
for the chemical potential at each of these extrema. Accordingly, the states corresponding 
to these extrema are self-consistent fragmented excited states of the condensate. 

The two maxima of E{ni) correspond to metastable states. Any variation of the boson 
occupation ni lowers the energy favoring a "decay" into either the GP ground state which is 
also the ground state of the present mean-field or to the state corresponding to the minimum 
of E{ni). The latter state is stable with respect to variations of rai, or briefly a stable state. A 
substantial energy barrier has to be overcome in order to lower the energy and to destabilize 
this state. Using equations (4-6) we can show that 

aril 

holds. This physically appealing relation adds insight into the understanding of the curves 
E/N depicted in fig.l. For large N the quantity /in — H22 is the energy gain obtained by 
removing a boson from the 0i-manifold and adding it to the 02-nianifold. /in — ^22 is thus 
the "driving force" for the flow of bosons between the two boson subsystems. At the extrema 
of -E as a function of the boson occupation one recovers (7). 

The orbitals 0i(a;) and 02 (a^) for the stable fragmented excited state discussed above are 



6 



depicted in fig.2A. For comparison also the orbitals (/? for the ground and for the first excited 
GP state are shown. The ground state GP orbital is as usual nodeless and consists of two 
humps related to the two wells of the potential. The first excited state GP orbital exhibits 
a node close to the top of the barrier. Each of the orbitals 0i and 02 of the fragmented 
stable excited state also possesses a node, but these nodes are located at distinct different 
sites. 4>i is rather localized at the deeper well and its node is close to the minimum of the 
other well. The reverse situation holds for 02- That these nodes are well separated can 
be understood by comparing the density per particles p — + 712102^)/-^ with the 

respective quantities pcp = for the GP states. These densities are shown in fig.2B 
together with the trap potential V{x). Since the bosons in our example repeal each other, 
the density would like to spread over space to reduce this repulsion. Indeed, in the ground 
state the density is substantial in between the wells. Due to the node of the excited GP 
state, the density vanishes close to the top of the barrier, and as a consequence the density 
is enhanced inside the wells. This inevitably leads to a substantial increase of the energy 
Ecp for the excited state. Fragmentation assigns different nodes to different orbitals and 
hence the density can penetrate the region of the barrier as seen in fig.2B. Consequently, 
the energy of the fragmented excited state is much lower than that of the GP excited state 
as seen in fig.l. 

It should be stressed that the above findings depend only weakly on the shape of the 
double-well potential. We have performed numerous calculations varying the coupling con- 
stant A, the bias A and other parameters of the trap potential. In particular, if we choose 
the potential to be symmetric (bias A = 0), the same conclusions as above can be drawn. 
Furthermore, enlarging the coupling constant A = AqA^ results in even more prominent 
minimum and maxima of the energy curves as a function of rii/N . 

Let us briefly discuss the short E/N curve in fig.l. We have obtained this curve by 
solving (6) starting from the GP excited state orbital, i.e., from n2 — N, and increasing rii 
continuously. As seen in the insert which shows this curve on an enlarged scale, the curve 
exhibits a maximum thus describing another fragmented metastable state. It is also seen 
that the curve has a boundary minimum at rii = 0, indicating that the corresponding GP 
excited state is stable. The energy barrier involved is, however, very low. We remark that 
at large values of A(> 8.0), this barrier disappears and we find a boundary maximum, i.e., 
the GP state becomes metastable. The E{ni) curve then acquires a minimum at ni 
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giving rise to a new stable fragmented state. 

In conclusion, self-consistent fragmented excited states of condensates exist. These states 
can be either metastable or stable with respect to the flow of bosons from one fragment 
(manifold of bosons residing in one orbital) to another one. In the class of examples studied 
here energies of fragmented states are much lower than those of the excited GP states. Due 
to the transparent physics behind populating several orbitals in excited states, we expect 
such states to exist also in two and three dimensions and in other trap potentials. 

The authors acknowledge useful discussions with Ofir Alon. 
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FIG. 1: Energy per particle E/N as a function of the fractional boson occupation ni/N for the 
coupling constant A = Aq = 2.5. For comparison, the energies Egp/N of the GP ground and 
first excited state are indicated. Two curves are shown, starting at the respective GP energies. 
The following parameters of the trap potential (8) have been used: A = 0.1, a = 0.04 and xq = 1.5. 
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FIG. 2: A: The orbitals (j)i{x) and (p2ix) corresponding to the stable fragmented excited state of 
fig.l. (for convenience {ni/N)^^'^(j)i(x) are shown) in comparison with the orbitals (/? corresponding 
to the GP ground and first excited state. 

B: The densities per particle p = (ni|(/>ip + n2\(p2\'^) / ^ of the stable fragmented excited state and 
PGP = Iv^P of the GP ground and first excited state. Also shown is the trap potential V{x) (for 
parameters, see caption of fig.l). The values of the potentials have been scaled by 1/20. 
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